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Abstract 

I calculate, at one loop in staggered chiral perturbation theory, the matrix elements of the 
complete set of five local operators that may contribute to B mixing both in the Standard Model 
and in beyond-the-Standard-Model theories. Lattice computations of these matrix elements by the 
Fermilab Lattice/MILC collaborations (and earlier by the HPQCD collaboration) convert a light 
staggered quark into a naive quark, and construct the relevant 4-quark operators as local products 
of two local bilinears, each involving the naive light quark and the heavy quark. This particular 
representation of the operators turns out to be important in the chiral calculation, and it results 
in the presence of "wrong-spin" operators, whose contributions however vanish in the continuum 
limit. If the matrix elements of all five operators are computed on the lattice, then no additional 
low energy constants are required to describe wrong-spin chiral effects. 

PACS numbers: 13.20.He, 12.39.Fe, 12.38.Gc 



I. INTRODUCTION 



The mixing of neutral B mesons provides a fertile area for precision tests of the Standard 
Model. The fact that the mixing is a second order weak process and is also suppressed by 
small CKM angles in the Standard Model makes it sensitive to new physics. In order to 
take full advantage of experimental measurements of the mixing, one needs to determine 
the hadronic matrix elements of the effective weak operators. For B mixing the relevant 
operators are local four-quark operators with Ab = 2, where b is 6-quark number, and the 
relevant states are B^ and B^ mesons or B^ and B^ mesons. A first-principle evaluation of 
such operator matrix elements is possible with lattice QCD. 

Lattice computations usually involve an extrapolation in light quark masses to the physical 
up and down masses, and always require an extrapolation in lattice spacing a to a = 0, the 
continuum. These extrapolations can be controlled by using a version of chiral perturbation 
theory that includes the effects of the discretization errors associated with the choice of 
lattice action. In two recent lattice calculations of B mixing [1, 2], staggered light quarks 
are combined with non-staggered heavy quarks using NRQCD [3] or the Fermilab action [4] , 
respectively. In such cases, the appropriate chiral theory is "rooted, heavy-meson staggered 
chiral perturbation theory" (rHMSXPT) [5]. 

In this paper, I calculate B mixing to one-loop order in rHMSXPT. Roughly speaking, I 
work to leading order in the heavy-quark expansion, although I do include the large I/uib 
effects: the hyperfine splitting of B and B* and the flavor splitting of Bs and Bd- This is a 
systematic approximation in the power counting introduced by Boyd and Grinstein [6] and 
discussed recently in Ref. [7] for the lattice calculation of heavy-light meson decay constants. 
If instead one prefers a power counting strictly in l/m^, which sets the splittings to zero, it 
is easy to take that limit of the results given in this paper. 

In the rHMSXPT calculation, it is important to take into account the exact form of the 
lattice operator used to approximate the continuum one. References [1, 2] construct the 
four-quark operators as the local product of two local bilinears,^ each formed from a heavy 
antiquark field and a light quark field with the "naive" lattice discretization. As proposed 
in Ref. [8] and discussed below, the naive field is constructed in turn from the simulated 
staggered fields (or more precisely, the naive propagator is constructed from the staggered 
propagator). Both the use of naive fields and the local nature of the four-quark operator 
influence the form of the corrections at one loop. 

It is not hard to understand the qualitative effect of the lattice locality of the four- 
quark operator. Because of lattice doubling symmetry, a (single-component) staggered quark 
field actually corresponds to the 16 continuum degrees of freedom of four "tastes" of four- 
component Dirac particles. On the lattice, the spin and taste degrees of freedom can be 
made explicit in position space by combining the 16 staggered components associated with 
an elementary hypercube [9]. For our four-quark operators, the two light staggered quarks are 
tied to the same space-time point, so that their spin and tastes are coupled. The coupling 
produces undesired contributions to the operator, with "wrong spin" and "wrong taste." 
These undesired contributions appear at 0{1) in the lattice spacing. Fortunately, in the 
matrix elements considered here, continuum SU{4) taste symmetry suppresses wrong-taste 
contributions and therefore wrong-spin contributions. On the lattice, SU (4) taste symmetry 



That is, all four fields are located at the same lattice point. 
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is violated at C(a^), so the undesired contributions come in at that order. Since C(a^) 
corrections appear at one loop in rHMSXPT, that is the order at which we find wrong-spin, 
wrong-taste contributions to B mixing. 

Similarly, it is clear that the effect of using naive quarks in the operators and interpolat- 
ing fields must be to sum over tastes, since the naive quarks have no explicit taste index. 
However, the details are non-trivial. It turns out that the heavy-light meson propagator 
is simply an average over the initial and final tastes, which are equal to each other. The 
three-point function involves a complicated sum over tastes of the staggered quarks in the 
interpolating fields and four-quark operator, and there is coupling between the spin matrices 
in the operator and the taste sum. These details play a key part in the discussion below. 

A calculation in rHMSXPT can be thought of as "staggering" the corresponding contin- 
uum calculation, which here would be in heavy-meson chiral perturbation theory. In fact, 
when the continuum calculation includes partial quenching effects, it is often possible to 
deduce the proper staggered version without having to recalculate explicitly any of the di- 
agrams (see, for example, Ref. [10]). In the current case, a partially quenched continuum 
calculation does exist [11]. However, the complications due to the naive-to-staggered trans- 
lation and the wrong spin-taste contributions make it necessary to perform the staggered 
calculation from scratch. Nevertheless, Ref. [11] is extremely useful here, and provides a 
check of the current results in the a — )■ limit. 

The B mixing matrix elements for any four-quark operator that can appear in the Stan- 
dard Model and in possible extensions such as supersymmetry can be written in terms of 
the matrix elements of the following five operators [12] 



where pairs of round or square parenthesis indicate how the color indices are to be con- 
tracted, and R and L are the right and left projectors: R = {1 + 75)/2 and L = (1 — 75)/2. 
Operators, Oi, O2 and C3 appear in the Standard Model, with Oi (which mixes with O2 
under renormalization) governing the mass differences of the neutral B eigenstates, AM^ 
and AMs- Operators O4, and appear in extensions of the Standard Model. Additional 
operators with R ■H- L can also contribute beyond the Standard Model, but parity implies 
that their mixing matrix elements in QCD are equal to those of the above operators. In 
addition to parity, Fierz transformations are needed in order to write the mixing matrix ele- 
ments of any four-quark operator with these quantum numbers in terms of those in Eq. (1); 
for a detailed explanation see Ref. [13]. 

I note that the corresponding projectors R,L in Ref. [11] do not have the factor of 1/2, 
so the operators there are differently normalized. Since in any case unknown low energy 
constants will enter in the chiral theory, this normalization difference is unimportant here. 

The fact that Eq. (1) is a complete set of operators for B mixing implies that wrong-spin 
contributions to the operators do not in fact lead to any new low energy constants in the 
chiral theory. Wrong-spin contributions to operator Oi merely lead to the appearance of the 
low energy constants associated with operators Oj^^i in the one-loop expression for the Oi 
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O2 

o. 



(hLq) [6Lg] 
ihLq] [bLq) 
(hLq) [bRq] 
(bLq] [bRq) , 
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matrix element. Thus, a staggered lattice calculation that computes the matrix elements of 
all the operators in Eq. (1) will not suffer from increased systematic or statistical errors due 
to the wrong-spin issue. Existing calculations [1, 2] study the matrix element Oi exclusively. 
In the case of Ref . [1] , the one-loop contributions of wrong-spin operators were not known at 
the time, so one presumably should include some additional systematic error in their result. 
In the case of Ref. [2], it was not possible to make a complete study of this effect because 
the matrix elements of the other operators were not computed. However, an associated 
systematic error was estimated. 

In the relevant staggered simulations, the fourth-root of the fermion determinant is taken 
in order to eliminate the four-fold multiplicity of tastes in the sea. The rooted theory then 
suffers from nonlocal violations of unitarity at non-zero lattice spacing [14, 15]. However, 
there are strong theoretical arguments [16-19], as well as other analytical and numerical 
evidence [20-24], that the local, unitary theory of QCD is recovered in the continuum limit. 
Furthermore, it is straightforward to take rooting into account in the chiral theory. One 
simply needs to multiply each sea quark loop by a factor of 1/4 [25, 26]. This can be done 
either by following the quark flow [27] to locate the loops, or — more systematically — by 
replicating the sea quarks rir times and taking = 1/4 in the result of the chiral calculation 
[17, 19]. Since I will need to work out the quark flows in any case, I use the former method 
below. 

The remainder of this paper is organized as follows. Section II briefly reviews the ba- 
sics of rHMSXPT, focusing in particular on those aspects that will be important here. In 
Sec. HI, I discuss the connection between naive and staggered quarks, and how it influences 
the structure of the four-quark operators and the interpolating heavy-light meson fields. 
The calculation of the one-loop diagrams is detailed in Sec. IV. I also briefly explain why 
taste-violations coming from mixing under renormalization do not need to be considered at 
this order. Section V compiles the final formulae for the chiral and continuum extrapolation 
of the matrix elements of the operators defined in Eq. (1); corresponding results for the B 
("bag") parameters are collected in Appendix A. I conclude in Sec. VI and make some addi- 
tional comments about existing and future lattice computations of B mixing. A preliminary 
account of the current calculation appears in Ref. [28]. 

Though I denote heavy quarks as b quarks and heavy mesons as B mesons throughout, 
the current calculation in rHMSXPT also applies to the local matrix elements in neutral D 
mixing, with the usual caveat that the omitted l/mg terms {Q is a generic heavy quark) 
are larger in that case. However, long distance contributions are presumably much more 
important in the D case [29] than in the B case. Such contributions are beyond the scope 
of this work, and are likely to be difficult to compute on the lattice. See however Ref. [30] 
for a lattice approach to long-distance effects in the kaon system. 

II. BASICS OF RHMSXPT 

Here, I give a brief summary of some of the basic features and definitions from staggered 
chiral perturbation theory, both of heavy-light mesons and of light mesons ("pions"). In this 
summary, I follow Ref. [10] fairly closely, but adapt the notation slightly, to make it more 
similar to that of Ref. [11]. The reader is referred to the literature [5, 10, 25, 26, 31, 32] for 
more details. For convenience in making connection to Ref. [11], I write the Lagrangian and 
do the perturbative calculations in Minkowski space. 



4 



Let Pg''"' be the field that annihilates the pseudoscalar meson containing a heavy quark b 
and a light quark q (the Bq for q = d), while P^^q^ does the same for the vector meson {B*q 
for q = d). To take advantage of heavy-quark spin symmetry, pseudoscalar and vector fields 
are combined in the field 



Hi") 



1 + 



7(Af) 



p*{b) I 



«75 



(A^) p{6) 



(2) 



which destroys a meson, while 



H 



' u p*(b)] , - AM) p{b)] 



(3) 



creates a meson. Here f is the meson velocity, and the (M ) on '^'^m) ^^"^ 75^^'' indicates that 

they are Minkowski-space matrices: '^^m) ~ 7°' 7/m) ~ '^7"'' ^^"-^ 75*^'' ~ 75' with and 75 
the Euclidean (Hermitian) Dirac matrices. The label q indicates the "fiavor-taste" index of 
the light quark in the meson. For n flavors of light quarks, q can take on 4n values. Later, 
I will write q as separate flavor (x) and taste (a) indices, q — )■ (x, a). 

Under SU{2) heavy-quark spin symmetry, the heavy-light field transforms as 



^ H 



(4) 



with S G SU{2), while under the SU^in)^ x SU^in)^ chiral symmetry. 



(b) 



(5) 



with U G SU (4n) defined below. We keep the light flavor and taste indices implicit here. 

The light mesons are combined in a Hermitian field For n staggered flavors, $ is a 

4n X 4n matrix given by: 



TT 

K- 



7r+ K+ 
S 



\ 



■J 



(6) 



I show the n = 3 portion of $ explicitly, and in fact detailed final results below will assume 
n = 3. Each entry in Eq. (A3) is a 4x4 matrix, written in terms of the 16 Hermitian basis 
elements of the Clifford taste algebra. It is convenient to take the generators of this algebra 
to be = 7* [9], where * denotes complex conjugation. Thus we write, for example. 



16 



(7) 



Th = {75, ^7^75, (Tf^uif^ < J^), 7m, 1} 
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with a^y = (i/2)[7^,7,,]. 

It is useful to divide the indices S into pairs of indices: H — )• (p, tp), where p labels 
the SO (4) representation (P,A,T,V,I) and tp labels the element within each representation. 
Thus tp runs from 1 to Np, where Np is the dimension of each representation (1,4,6,4,1, 
respectively). We then can write 



No 



(9) 



The component fields of the flavor- neutral elements of $ (namely Up^tp, Dp^tp, 
real; the other (flavor-charged) fields (vr^t^, K^^^, . . . ) are complex. 
The mass matrix is the An x 4n matrix 



are 



M 



(mj 

rudl ■■■ 

rrisl ■■■ 

V ; ; ; 



(10) 



where the portion shown is again for the = 3 case. 

From $ one constructs the unitary chiral field E = exp[z<l>//], with / the tree-level pion 
decay constant. In our normalization, / ~ = 131 MeV. Terms involving the heavy-lights 
are conveniently written using cr = a/E = exp[i$/2/]. These fields transform trivially under 
the SU{2) spin symmetry, while under SU{An)L x SU^in)^ we have 



S - 



(11) 
(12) 



with global transformations L G SU (4?^)^^ and R G SU{An)R. The transformation U, defined 
by Eq. (12), is a function of $ and therefore of the coordinates. 

It is convenient to define objects involving the a field that transform only with U and U^. 
The two possibilities with a single derivative are 



[a'^dpa - adpa'^] . 



(13) 
(14) 



transforms like a vector field under the SU{An)L x SU{An)]i chiral symmetry and, when 
combined with the derivative, can form a covariant derivative acting on the heavy-light field 
or its conjugate: 



iH('%), = ^ Spiff + 



0pH% = l^fH^} ^ 5pFf - ^YfH^} , (15) 
with implicit sums over repeated indices. The covariant derivatives and Ap transform under 
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the chiral symmetry as 

^ UA^,U"^ . (16) 
We can write the leading order (LO) chiral Lagrangian as 

^LO = -Cpion + ^UL , (17) 

where £pion is the standard staggered chiral perturbation theory (SXPT) Lagrangian for the 
light-light mesons, and £hl is the contribution of the heavy-lights. In Minkowski space, we 
have 

Cpion = ^Tr(a^Sa'^St) + i/i/2Tr(A^S + A^St) 

-^(f// + ^/ + 5/ + ...)'-«'V, (18) 
-V = CiTr(e5Se5St) + C3^ J][Tr(e.Se.S) + /i.c.] 

+ ^"4^ Yj^^iU^lliUr^T.) + h.c\ + Ce 5^ Tr(e^,Se^,St) 

+ C2y^ $^[Tr(e.S) Tr(e.S) + h.c\ + Csa^ 5^[Tr(2e.e5S) Tr(ze.e5S) + h.c\ 

V V 

+ C5y^$^Tr(e,S)Tr(e.St) + C5A^ J]Tr(ze.e5S)Tr(z^.6S^) , (19) 

V V 

£hl = -^ Tr(F(%(^) t;-&) + Tt(h'^''^ H('%,^^i'''^A,) . (20) 

Here Tr denotes a trace over flavor-taste indices and, where relevant, Dirac indices. The 
product H^'''^ H^^^^ is treated as a matrix in fiavor-taste space: {h'^'^^ H''^^)qqi = H^^^ H^'l'\ The 

covariant derivative & acts only on the field immediately preceding it. For convenience, I 
work with diagonal fields {U, D, . . . ) and leave the anomaly (/«o) term explicit in Eq. (18). 
We can take rriQ — )■ oo and go to the physical basis (vr", rj, . . .) at the end of the calculation 
[33]. 

At tree level, the light-light meson composed of quarks of flavor x and y, and with 50(4) 
taste representation p, is 

^i,p = /^K- + my) + a'Ap . (21) 

Here Ap is the taste splitting, which can be expressed in terms of Ci, C3, C4 and in 
Eq. (19) [25]. The residual 5*0(4) taste symmetry [31] at this order implies that the mesons 
within a given taste representation are degenerate in mass. 

I now list some key expressions from the Feynman rules given in Ref. [10], but adapted 
to the current notation. Using separate indices for flavor {x,y) and taste {a,a',c,c'), the 
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(quark-line) connected pion propagator in Minkowski space is 



J conn 



a cc' 



tp 



(22) 



Here, I have used the fact that V'^p is Hermitian is to replace the complex conjugation in 
Eq. (9) by interchange of indices on the right-hand side, which will be convenient later. 
Similarly, the disconnected (hairpin) propagator is 



(23) 



where 



xx,yy 



Ml) 



' (p2 - + le) (p2 _ + le) {p^ - M^^ + te) {p^ - M^, ^ + ze) 
with the hairpin strength b'p given by 



(24) 



51 



4m2/3, 
0, 



p = V (taste vector); 

p = A (taste axial- vector) ; 

p = I (taste singlet); 

p = T,P (taste tensor or pseudoscalar) 



(25) 



X and Y denote valence mesons made from xx or yy quarks, respectively, with Mx and My 
their masses. For the sea mesons, the masses Mjj and Ms do not include the mixing effects 
of the hairpins. The re-diagonalized states after including the hairpins are t] and t]'. For 
concreteness I have assumed the 2+1 case: m„ = ma. 
The propagators for the heavy-light mesons are 



p*(h) p*{b)]\/i^\ 



2{vk + ie) ' 
2(vk - A* +ie) 



(26) 
(27) 



where A* is the B*-B mass splitting. The B B* n vertex (including the i from exp{iC)) is: 

(28) 



^ (pm p*{h) _ p*{b)\ p(h)\ 

J \ xa ^yc fixa yc J 



yx 
ca ' 



where the repeated indices a, x, y, c, and /i are summed. 

For B mixing, we also need corresponding fields that destroy and create mesons with b 
quarks, i.e., i?o-like and i?Q-like mesons. These fields and their interactions can be obtained 
from the previous ones using charge conjugation [34]. The individual meson fields are indi- 
cated by, for example, P^^q'' and Pq''\ (Note that the light quark label q does not distinguish 
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between quarks and antiquarks.) The combined fields are 



p*{b) _|_ 
'{My IJ.,q ^ 



{M)p{b) 



H 



(b) 



' u p*(fe)t , - AM) p(b)t' 



(29) 
(30) 



The propagators for the P^^g^ and Pq"' fields are the same as those for the P^'^q' and Pg 



*ib) 



,ib) 



fields, Eqs. (26) and (27). The BB*tt vertex is 

f pib) p*(fe)t _ p*(b) p(fe)t A fjf^^y^ 
\ xa fiyc fixa yc j ^ ca ' 



(h 
f 



(31) 



III. TRANSLATING FROM NAIVE TO STAGGERED QUARKS 

The naive light quark action may be rewritten as four copies of the staggered action: 

^{x) = Q{x) x{x) ; fi(x) = 7o° 7r 72^ 73^ (32) 

where \l/(x) is the naive quark field and x{^) is a "copied" staggered field, with each Dirac 
component Xi separately having the staggered action. I call "copy symmetry" the SU (4) that 
acts on the copy index i. Unlike the SU{4) vector taste symmetry, which acts on individual 
staggered fields (written in the spin-taste basis) and is exact for an interacting theory only 
in the continuum limit, copy symmetry is an exact lattice symmetry. Thus the propagator 
of a copied staggered field is 

{xi{x) Xi'iy)) = {x{x) x{y)), (33) 

where x is the normal (uncopied) staggered field. This implies 

{^ix)^y)) = nix)n\y){x{x)xiy)) . (34) 

In the simulations using staggered quarks, the naive field is never constructed per se; instead 
Eq. (34) is used to translate staggered propagators into naive propagators [8]. 
An interpolating field 'H(x) for a Bg meson is 

n{x) = b{x) 75 *(a;) = b{x) 75 fi(x) x{x) ■ (35) 

I assume that in practical applications 'H(x) will always be summed over a time-slice, either 
explicitly, or implicitly by using translation invariance. 

To leading order in a, b{x) varies smoothly (up to gauge transformation) between neigh- 
boring spatial sites, but x does not, due to taste doubling. On the other hand, in the 
spin-taste basis, which we arrive at by summing the staggered fields over hypercubes, the 
staggered fields are smooth on the doubled lattice. We are thus led to focus on the average 
of ^{{x) over a spatial cube. Let x = {t,x) with x = 2y even, and let r] = {r]o,r]) be a 
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4- vector with all components or 1. For t even {t = 2r) the averaged field is 
n^^-'Kt^x) = ^5^%aj + r7)75fi(2r,r7)x(2r,22/ + 77) 

= ^6(t,a;)75^fi(T7)x(2r,2t/ + r7) 

= ^ bit, x) 75 J2 [^(^) ^(2^ + ^0' 2?/ + ^) + 

V 

+ (-!)"« n{r]) x(2r + r^o, 2?/ + ij)] . (36) 

Inserted gauge links for gauge invariance of point-split quantities are implicit. For t odd 
(t = 2r + 1), the result is the same except the term on the last line of Eq. (36) changes sign. 
Using the fact that (—1)''" fl{ri) = 7570 ^{rj) 7075, it is not hard to see that this second term 
just gives the usual staggered oscillating (in time) state with opposite parity. I have dropped 
higher order terms in a coming from the variation of the heavy-quark field over the cube. 
For simplicity, we simply assume from now on that all components of x are even {x^ = 2y^), 
and that the oscillating state has been removed by the fitting procedure. We then have 

n^'^'^Hx) -> ^&(x)75 Y.^{v)x{2y + V) ■ (37) 

V 

We now convert to a spin-taste basis for the staggered fields. The standard construction 
for a single staggered field is [9] 

= ^ 5^r]"'^(r/)x(2y + r7) , (38) 

where a is a spin index and a is a taste index. As is well known [35-37], this decomposition 
is correct only to lowest order in a and generates a spurious 0{a) term in the spin-taste 
action, but it is good enough for our purposes. Here we need a copied version: 

Qriy) = l 5^^^"'^(r7)x.(22/ + r/) • (39) 

With spin indices implicit, Eq. (37) then becomes 

■H^''''\x)^h{x)^,qt{y)6t , (40) 

where repeated indices are summed. With Eq. (33), this implies that the contraction of H 
with "H^ [i.e., the heavy-light propagator) is automatically averaged over tastes: 

{n{x) H\x')) ~ \ {Kx)^,q\y) riy'HKx')) , (41) 
where a sum over a is implicit. 
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Analysis of the four-quark operator is more complicated because the two bilinears from 
which it is constructed are not separately summed over space; only the four-quark operator 
is summed. However, we can write it in terms of separately summed bilinears by using the 
identities 

^J2^T{n{v)Kn^{v)K) ii[n{^')Kn\^')K) = <5,,,, , (42) 

^tT{Q{r])KQ\r])K) Q{r]) = KQ{ri)K . (43) 

Here K is any of the 16 independent Hermitian matrices in Eq. (8), which obey = I. 
We get, for operator On = &r„^ &r^^: 

(^i^'Ha^) = ^^6(t,2|/ + r/)r„v]>(t,22/ + r7)6(t,2y + r/)r;vl/(t,22/ + r7) 

V 

-> ^ J2(bTnKql hV'^Kqi) K,,K,e • (44) 

K 

Here we have dropped the "wrong parity" part, which does not contribute if oscillating terms 
are removed by the fitting procedure. Note that contributions with K ^ I have incorrect 
spin {TnK (g) T'^K instead of r„ ® r'„), and coupling of taste (c, d) and copy {k, i) indices. 

There are 0{a) and higher corrections to Eqs. (40) and (44), coming from the variations of 
the heavy quark field over the spatial cube and from corrections to the spin-taste construction 
in position space. As discussed in Sec. IV C, however, such terms do not contribute to non- 
analytic terms in chiral perturbation theory until NNLO. 

Using Eqs. (40) and (44), copy symmetry (Eq. (33)) implies 

^^(a^)t Qiav) ^(a^)t^ ^ {D^^ D'^^) K^aKde , (45) 

where D'^'^ is the quark propagator (in a given background) for taste a into taste c. If taste 
symmetry is exact, {D'^°-D'^'^) oc Sacked, and only the correct spin {K = I) contributes. Thus 
the desired matrix element will be obtained in the continuum limit. 

At one loop, however, taste violations allow (D'^^D^'^) not to be proportional to 6acSed, 
and terms with incorrect spins can contribute. For example, the taste-violating hairpin with 
vector taste can give a term proportional to ica^de — laded- Then, since tr^jf^K) = 4:SK,'y^, 
the spin of the operator is r„7^ ® r^7^ instead of r„ ® r'„ 

Equation (44) may be simplified by taking advantage of the exact S'0(4) taste symmetry 
of the staggered chiral theory at one loop. Within any S'0(4) multiplet k G {P, A,T,V, 1} 
with dimension iV^ > 1 {e.g., the vector-taste multiplet V with Ny = 4), the value of all one- 
loop diagrams would be unchanged if we replaced any multiplet element in the taste factors 
KckKdt with another element from the same multiplet.^ We therefore write K = F'^'*", where 
tfj_ labels the element within multiplet k. Replacing the sum over K in Eq. (44) with a double 



^ Although indices k and £ are copy, not taste, indices at this point, they will become taste indices d la 
Eq. (45) shortly. 
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sum over k, t^, and using the 5*0 (4) symmetry, we then have 

Within a given multiplet, this decouples the sum over spins from the sum over tastes. 

Thus, for a given continuum operator 0„, the plan is to calculate the one-loop diagrams 
for each of the operators 

o:;^5^6r„r'''*''=g^6r;r«'*«g,^ (47) 

between external (interpolating) fields ("H^)" and {T-OYj, where, from Eq. (40), 

= 9- 75& • (48) 

Each diagram for operator O'^ then gets an additional factor coming from Eqs. (46) and 
(40), where 

^'^^I^ErS^rf'^^.A,. (49) 

Whether explicitly indicated or not, all repeated indices will then be summed; this includes 
taste and copy indices {a,c,d,e and i ^ J ^ £^ cis well as the indices with dual, spin-taste 
meaning {n.tf^.t'^. 



IV. CALCULATION OF ONE-LOOP DIAGRAMS FOR 0„ 
A. Procedure 

We now set up one-loop rHMSXPT for the operators described above. We follow Ref. [11] 
as much as possible, but must take into account the complications of copy and taste indices. 
It is convenient first to express the operators (9^, given in Eq. (47) in terms of the basis of 
Eq. (1). From the relations among operators listed for example in Ref. [13] we find 
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Of = Ol 
O^ = -8O2 - 8O3 , 

of = -QOi , (50) 
C]" = 8O2 + 8O3 , 
Oi = Ol, 

Ol = O2 , 

O^ = -Ol , 

Ol = - 4O3 , (51) 

Ol = Ol , 

Oi = O2 , 
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ox 



-Ox , 

-A02 - 2C3 
03, 



(52) 



Of 

ol 
o\ 

o{ 



-04 , 

, 

-205 , 

, 



(53) 



Ol 

ot 
ol 
o\ 
o{ 



-o, , 

-204 , 

, 

-204 , 
. 



(54) 



The chiral representatives of the standard operators on the right hand side of Eqs. (50) 
through (54) are given in Ref. [11]. There, the only relevant quantum number of the light 
quarks is their flavor, and both bilinears have the same flavor, which is labeled q. Here we 
also need to label the taste and, for the moment, the copy index of the light quarks, and these 
are not in general the same for both light quarks in the operator. So we adopt the notation 
g — )■ X, c, fc, where x labels the quark flavor only, c (or other letters near the beginning of the 
alphabet) labels the quark taste, and k (or other letters near the middle of the alphabet) 
labels the quark copy. From [11], we then have 



01 



xck:xd£ 



O: 
O 



/3i 



xck;xd£ 
2(3) 

xck;xd£ 
4(5) 



x,c,k \ / x,d,£ \ / x,c,k V 



xAl 



rpm\ 



/32(3)K^^"") , 

/ x,c,k V / x,d, 



rp(b) 



/34(5) 

2 

I ^4(5) 



x,c,k \ J x,d,£ 



/ x,c,k 



x,c,k \ J x,d,l 



crt p*(b)\ 



x,c,k 



^p*{b),ii. 



x,d,l 



(55) 



The method used to obtain these operators in Ref. [11] is a standard spurion analysis. The 
factors of a and o"^ are present in order make the light-quark spurious, which transform 
by left or right chiral rotations in Eq. (1), into objects that transform with U (defined in 



Eq. (12)) and can combine with the heavy-meson fields H^^^ and Hg"' to make invariants. 
Although many insertions of Dirac matrices are possible in forming the invariants, they all 
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reduce down to the simple forms in Eq. (55) when expressed in terms of P'^^^''', p(*\ Pl!'^^\ 
and P^^^\ As pointed out in Ref. [34], this follows from heavy-quark spin symmetry, which 
relates the amplitude for B-B mixing to that of B*-B . 

In Eq. (55), I have used the fact that we are only interested in the parity-even part of 
these operators to set the two coefficients Detmold and Lin call (3^^) and /34(5) simply to 
/34(5)/2. So where they have /34(5) + /34(5) we will have simply /34(5), and similarly for (3'^^^^ 
and /34(5). We use the Latin O for the chiral operators to distinguish them from quark- level 
operators O. The external interpolating B and B fields are taken to be, respectively, 

P^'L and PiSi-. (56) 

Strictly speaking, there should be a second set of terms on the right hand sides in Eq. (55) 
in which the pairs of indices c, k and d, I are interchanged. These come about because the 
light quark field in either of the bilinears of the operator can be the one that creates the 
light antiquark in the B or annihilates the light quark in the B. However, since Eq. (49) 
is symmetric under this interchange, the extra terms give identical results to the ones we 
have already, and therefore can be dropped at this point. Note that the fii are low energy 
constants with arbitrary normalization. 

We can now use copy symmetry to simplify the equations, and ultimately eliminate the 
copy indices entirely. However, this cannot be done without taking into account the quark 
fiow through the diagrams: copy symmetry works at the level of the light quark propagators, 
and not on meson propagators 'per se. The point is that a given light quark in an external 
meson field can end up in either a (6)- or (6)-labeled meson field in Eq. (55), and the copy 
symmetry has a different effect on the diagram in the two cases. If the light quark in the 
external P^*] ^ field contracts with the quark in a (6)-labeled meson field (which has taste 
and copy indices c, k), then the combination of copy symmetry (which gives 6ik) and the 6ai 
in Eq. (49) forces k = a. On the other hand, if the same light quark contracts with the 
quark in a (6)-labeled meson field (indices d,i), then we have i = a. So P^ will end up 
with a different taste structure in the two cases. For convenience, I prefer to rename the 
taste indices (c -H- d) in the second case so that P^ remains the same. In doing so I adopt 
the convention that the taste-a quark in the external (6)-labeled field always contracts with 
the taste-c quark in the operator, but that taste-c quark may be in either the (6)- or the 
(6)-labeled meson field of the operator. Similarly, the taste-e quark in the external (6) field 
always contracts with the taste-ci quark in the operator. 

Keeping the above points in mind, we now eliminate the copy indices completely from 
the calculation. With the exception of the wave-function renormalization diagrams, which 
of course do not involve the four-quark operators at all, the final procedure is as follows: 

1. For a particular operator On of interest, we first write the related operators O'^ as 
linear combinations of the standard operators, following in Eqs. (50) through (54). 
It will be convenient then to define and (3n'^^ as the /3 and (3' corresponding to 
operator O'^. Table I for /SiC"* follows immediately from Eqs. (50) through (54). 

2. We then calculate the chiral diagrams, using copy-free versions of Eq. (55) for the 
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TABLE I: Values of (3^ ■ For pT\ simply put primes on all entries in the table, with the under- 
standing that f3'i = Pi (see Eq. (55)). 



n 




P 


A 


T 


\/ 


/ 


1 


Pi 


-8/32 - m 


-6/3i 


8/32 + 8/33 


/3i 


2 




-/3i 


-2/32 - 4/33 


/3i 


/32 


3 




-/3i 


-4/32 - 2/33 


/3i 


/33 


4 


-/34 


-2/35 





-2/35 


/34 


5 


-/35 


-2/34 





-2/34 


/35 



chiral operators, namely 



0xc;xd 

^xc;xd 
2(3) 



/3i 



a 



p(b) 



x,d 



a 



a 



p*ib),fi 



xA 



[or c ^ c/] , 



O 
O 



V / x,c \ / \ / x,c \ / x,d 



4(5) - 2 
/^4(5) 



cr 



+ 



/ x,c V / x.d V / 2:,c V / X, 



The external interpolating fields are 

Pil) and Pil)t 



(57) 
[or c ^ c?]. 

(58) 



For a diagram with a given quark flow, one should use either the explicit version of 
the operators in Eq. (57) or the alternative c ^ d forms to ensure that the external 
taste-a light quark contracts with the taste-c light quark in the operators (which also 
guarantees that tastes e and d contract). 



3. Each diagram is then multiplied by the overall factor 

P, 



16N, 



/ ^ ca de ' 



(59) 



and the repeated taste (a, c, d, e) and spin (k) indices are summed. 



B. Calculation 

We are now ready to compute the one-loop diagrams in rHMSXPT using the Feynman 
rules given above in Eqs. (22) through (28) and Eq. (31). Because of the complications 
due to taste and copy indices, it is not possible in general simply to modify the continuum 
results of Ref. [11] to insert staggered corrections, as in Ref. [10]. We must calculate most 
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FIG. 1: Meson- level tadpole graphs for the B-B mixing matrix element. In (a), the pion fields 
that are contracted both come from the same factor of a in Eq. (57), while in (b) they come from 
different a factors. For definiteness, we take the right external line in each diagram to be the 
incoming B meson, and the left external line, the outgoing B. This means the left line is the 
contraction with the P^^')'^ field in Eq. (57), while the right line is the contraction with the P^^^ 
field. Diagram (a) is "factorizable" into a product of the left- and right-hand parts (the right-hand 
part is trivial here). Diagram (b) is "non-factorizable." There is another diagram equivalent to (a) 
in which both pions come from the a associated with the right-hand line. 



chiral diagrams from scratch. The exception is the wave-function renormalization (param- 
eterized below by the function W), which is simple enough that the modification process 
("staggering") works. In the wave-function case, the naive-to-staggered translation gives no 
complications because, from Eq. (41), it only requires setting initial and final tastes equal 
and averaging over them. The taste-averaging has no effect because discrete taste symmetry 
(shift symmetry) implies that a two-point function is in any case proportional to the identity 
in taste space. ^ 

In addition to wave-function renormalization, there are two types of one-loop diagrams: 
tadpole graphs. Fig. 1, and sunset diagrams. Fig. 2. These are parameterized by functions T 
and Q respectively. Contributions from incorrect spins can enter in the tadpoles and sunset 
diagrams; we call such contributions T, and Q. The complete matrix elements are given by 

(BllOtlBl) = A (l + ^^k±]^ + ri^) + ri^) + Q« + Qi^)) + analytic terms. (60) 

and 

(S^IO,.") . (l . , r.". . r<»') . K (Q<"> . Si"') . a„a„.ic .... 

(61) 

for n = 2,3,4,5. Nonrelativistic normalization, which is standard in heavy-light chiral 
perturbation theory, is assumed for the states {B^\ and \B^) in these expressions. With 



^ With the exact momentum space taste construction [35], this statement is true to all orders in a, as can 
be seen most easily by using the formulation of shift symmetry in Ref. [19]. The construction is built into 
staggered chiral theory, so the statement is also true to all orders in rHMSXPT. The fact that we have 
used the (approximate) position space taste construction [9] in the translation from naive to staggered 
operators is irrelevant, since for our purposes only the lowest order translation is needed. 
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FIG. 2: The meson-level sunset graph for the B-B mixing matrix element. Conventions are as in 
Fig. 1. 



relativistic normalization, an extra factor of Mb^, the mass of the Bx meson, would appear 
on the right-hand sides. 

In the partially quenched 2+1 (m^ = 7^ m^) case, staggering the result for Wgi in 
Ref. [11] gives: 



xb 



bx 



P 



16 



Xil,p 



+ 



4?({<'.};{fv})||^ 



(62) 



where the index p runs over the taste representations (P,A,T,V,I) with degeneracies Np, and 
§ runs over the sea mesons m, c?, s. The function % is equivalent to the integral H{m, A) 
defined in Eqs. (A3) and (A12) of Ref. [11]. The subscripts on H implicitly give the meson 
mass m in H[m, A) by specifying its flavor and taste. The flavor is indicated either by giving 
the flavor of the two quarks in the meson, as in x§ in the flrst term in Eq. (62), or by giving 
the name of the meson, as in the remaining terms, where X refers to the meson made of two 
light valence quarks xx {mx is its mass). The superscript on l-i is the second argument of 
the function H{m,A). It is the mass splitting between the heavy-light vector meson in the 
chiral loop and the external heavy-light pseudoscalar meson. In addition to the hyperflne 
splitting A* = Mb* — Mb, it includes a light flavor splitting whenever the light flavor of 
the vector meson in the loop is different from the the external flavor. For the flrst term in 
Eq. (62), the vector meson has flavor S so the splitting is 6sx = Mb^ —Mb^ = 2Ai/i(ms — m^), 
where Ai and /z are low energy constants. The constant Ai comes from heavy quark effective 
theory, and /i relates light meson masses to quark masses, Eq. (21). 

For comparison, the function % is the same (up to constants) as the function J introduced 
in Ref. [7] (Eq. (6.17)). Similarly the function Xj p deflned in [11] and used below is the same 
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up to constants as the function i^m'jp) used in Refs. [5, 10, 26]. The relations are: 



^^t = -y|^'^K-,p'^)' (63) 

^^.-p = I^^(^L) = Yl^2<pHml/Al) . (64) 
In the hmit of no sphttings, 

<P = -3^^.,P = -ll^^(^L) • (65) 

If one wants the B mixing result in the strict l/ms power counting in which the splittings 
are set to zero, one can simply use Eq. (65) for Ti everywhere below. 

The (Euchdean) residue functions i?'"'*^^ and Z}^"''^^ in Eq. (62) are defined by [26] 



<'"(W,M) 



J 3,1 



Dt''\{mh {/^}) - -^<''^({^}' {/^}) • (66) 



The mass combinations appearing as arguments of these functions in the 2+1 partially 
quenched theory are 

{fi} = {mi,mh} . (67) 

The tastes of these mesons (J, V, or A) are indicated explicitly in Eq. (62). 

The staggered heavy-light wave function renormalization is also calculated in Refs. [5, 10] 
for the case where the heavy-meson splittings are neglected; the result of adding in those 
splittings as explained in Ref. [7] agrees with Eq. (62). 

The tadpole and sunset contributions are more complicated, and we must follow the 
procedure outlined at the end of Sec. IV A. All four taste indices (two from the interpolating 
fields and two from the light quarks in the four-quark operator) enter in a non-trivial way, 
and shift symmetry does not require that all be equal. Indeed, taste symmetry violations 
arising ultimately from high-momentum gluon exchange can in general make one pair of 
tastes indices different from the other pair. However, if there are parts of a diagram that 
give simply a tree-level heavy-light propagator, the average over taste implied by Eq. (41) can 
suppress the taste-changing interactions. On the other hand, in some diagrams for wrong- 
spin contributions, the overall factor, Eq. (59), can project onto particular taste- violating 
internal pion propagators. The bottom line is that one must calculate the tadpole and sunset 
diagrams from first principles, and not simply try to stagger the continuum result.^ 



^ The rules from Rcf. [10] for staggering a continuum result would apply unchanged if, for example, we had 
set the two tastes in the four-quark operator equal and averaged over them, while either fixing the external 
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FIG. 3: Connected quark flow tadpole diagram. This is a contribution to Fig. 1(a). The thin 
lines denote staggered light quarks, and the think, gray lines denote heavy quarks. Indices shown 
(a, /, c, d, e) are taste indices. The left pair of touching filled black circles represents the light and 
heavy quark in the P^^^tj right pair represents the light and heavy quark in the P^''). Taste 
conservation for the right-hand B propagator forces d = e. 

I start with the tadpoles. I call a diagram "connected" or "disconnected" based on 
whether the internal pion propagator is connected or disconnected in the quark flow sense 
of Eqs. (22) and (23) above. The only connected quark flow possible for Fig. 1 (a) is shown 
in Fig. 3. This is a diagram with a sea quark loop; the sea quark has taste /, which is 
summed over. Taste conservation for the right-hand heavy-light propagator (the p(^)p('')t 
propagator) gives a factor of 6de- When combined with Eq. (59), this implies that F^,*^ = /, 
so this diagram has no wrong-spin contributions. The taste factor from Eq. (22) is simply 
TafTfc = Sac, which combines with the 6ca from Eq. (59) (using F^^f^ = J) to give a factor of 
4. The sum over tp then gives Np, the degeneracy of representation p. Including the factor 
of 2 from the equivalent diagram with the loop on the right side of Fig. 3, and the factor of 
1/4 for the rooted sea-quark loop, we get 

where S runs over the sea quarks m, (i, s. Note that the result is independent of n. Since 
wrong-spin contributions to this diagram are absent, and (5n'^ = fin for all n. Fig. 3 is simply 
proportional to /3„, which is is factored out of T in Eqs. (60) and (61). The difference in chiral 
structure between the left-right operators 04,5 and the left-left operators Oi,2,3 in Eq. (57) 
is not relevant because the even terms in the expansion of a and cr^ are the same, and both 
TT flelds in Fig. 1(a) come from the same a or factor. 

There is also a disconnected contribution to Fig. 1(a). The quark flow diagram is shown 
in Fig. 4. Again the 5de from the right-hand B propagator means that only the correct spin 
contributes. The taste factor from Eq. (23) is VafV = ^ac, which gives a factor 4 after using 
Eq. (59) and a factor of Np from the sum over tp. In this case only the I, V, and A channels 
contribute; see Eq. (25). In the singlet case, we can use the fact that M!^, j ^ ml for large 
mo, and take the limit mo — ?• 00, resulting in one less pole in the denominator of Eq. (24) and 
an overall factor of 4/3. In the vector and axial channels, we simply get a factor of 4 from 
Np. This gives the standard ratio 1/3 between the strength of the singlet and vector or axial 



tastes to any one value, or averaging over them. Such a definition of the operator could be arranged if we 
constructed point-split bilinears with the desired tastes within a hypercube, and then multiplied two of 
them appropriately. In that case, there would not be any wrong-spin contributions. 
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FIG. 4: Disconnected quark flow tadpole diagram. This is a contribution to Fig. 1(a). Taste 
conservation for the right-hand B propagator forces d = e. 



hairpins, as seen in [10]. One may expect this standard ratio in any meson diagram, such as 
Fig. 1(a), that is unaffected by the comphcations from wrong spins or the naive-to-staggered 
translation. The result is 



' x,Fig.4 




I 2 c-/ 

+a dy 



dixy 



dmlr 



E 4S({mS};W})x,, 



+ {V^A) 



(69) 



Again, since wrong-spin contributions to this diagram are absent, the result is independent 
of n. 

We now turn to contributions to the non-factorizable diagram. Fig. 1(b). The connected 
contribution is shown in Fig. 5. This is the first diagram in which the quark flow connects 
the taste-a quark in the external Pifi field with the light quark in the P*-^-* field of the 
operator. According to the discussion above, this requires that we use the c <H- c? versions of 
the operators in Eq. (57). The combination of taste matrices F^^'' Tec'' from the connected 
propagator Eq. (22), and T'^^'^T'^'^" from the overall factor, Eq. (59), appears in several 
diagrams. It is therefore useful to define 

C{k, p) = Yl (r'^'*'=F^'*''F^'*'= F^'*") . (70) 

The factors of A^^^ and Np have been included in the denominator for later convenience. 
Completeness of the 16 matrices F*^'*" implies that ( satisfies the normalization condition, 

5^Ar,C(«:,p) = 16 5,,, . (71) 



Values of ({k, p) are given in Table II. 
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FIG. 5: Connected quark flow tadpole diagram. This is a contribution to Fig. 1(b). Taste-violations 
on the connected pion line allow for a ^ c and d ^ e. 



TABLE II: Matrix of C(k, p) values. 





P 


A 


T 


V 


/ 


p 


1 


-1 


1 


-1 


1 


A 


-1 


-1/2 





1/2 


1 


T 


1 





-1/3 





1 


V 


-1 


1/2 





-1/2 


1 


I 


1 


1 


1 


1 


1 



In terms of ({k,,p), the correct-spin and wrong-spin contributions of Fig. 5 are then 

(72) 




ri:pl5 = E ( E C(«:, P) X.J , (73) 

where the upper sign is for n = 1, 2, 3 and the lower sign is for n = 4,5. The difference in 
chiral structure of Eq. (57) between O4 5 and Of 2 3 in Eq. (57) matters for Fig. 1(b), because 

the two pion fields come from different a or factors. The correct-spin contribution 7^^p|g 5 
comes from the k = I term in the sum, while the k ^ I terms give the wrong-spin piece 
T^^pjgg. The quantities fi^^ are listed in Table I. 

The final tadpole diagram is the disconnected contribution to Fig. 1(b), shown in Fig. 6. 
In this case, the taste structure of the disconnected propagator Eq. (23) is FacTg^*'', which, 
when combined with the taste matrices in the overall factor, Eq. (59), gives tr^(r'^'*''r'''*'') = 
16 5i^p5t^tp- The hairpin propagator is non-zero only in the /, V, and A channels. So 
we get a singlet contribution to the correct-spin operator = p = I) and taste-violating 
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FIG. 6: Disconnected quark flow tadpole diagram. This is a contribution to Fig. 1(b). Taste- 
violating hairpins allow a ^ c and d ^ e. 



contributions for two of the wrong-spin operators {k = p = V^ A). The result is 

dXxj 



x,Fig.6 



-fin) 
' x,Fig.6 



3/2 



(74) 



dX 



xy 



-E 



(3)- 



■j,v 



(75) 



dl 



X,A 



Again the upper sign is for n = 1, 2, 3 and the lower sign is for n = 4, 5; the reasoning is the 
same as in Eqs. (72) and (73). 

The sunset diagram, Fig. 2, is very similar to the tadpole contribution that connects 
the incoming B and outgoing 5, Fig. 1(b). Again, there are two quark flows: Fig. 7, a 
connected graph similar to Fig. 5, and Fig. 8, a disconnected graph similar to Fig. 6. The 
taste structures of the sunset graphs are identical to the corresponding tadpole graphs. As 
in Fig. 5, the c ^ d version of Eq. (57) is used in Fig. 7. The main difference between 
the sunset and tadpole graphs is the actual integral, which here involves two heavy-light 
propagators and factors of gs'Sn, and so gives g'^,^^'H^* instead of X. In addition, the a 
and (j^ matrices in Eq. (57) are all set to 1 in the sunset case, so there is no difference in 
overall sign between the results for operators 1, 2, 3 and those for operators 4, 5. Otherwise, 
everything is the same as for the tadpole case. I find, for Fig. 7, 



*^x,Fig.7 



Qin) 



'^9b*B-k 

■ 16/2 



%Ef|^E^.c(K,p)7^^; 



16/2 



(76) 
(77) 
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FIG. 7: Connected quark flow sunset graph corresponding to Fig. 2. Taste- violations on the 
connected pion line allow for a ^ c and d ^ e. 




FIG. 8: Disconnected (hairpin) quark flow sunset graphs corresponding to Fig. 2. Taste- violating 
hairpins allow for a ^ c and e. 



The disconnected sunset graph, Fig. 8, gives 



Q 



in) 

x,Fig.i 



Qin) 
^^x.Fig.S 



3/2 



E 4S({<};{M)^^^; , (78) 



ie{Mf)} 



dHxy 



(3)- 



(79) 



X,A 



C. Other Possible Taste-breaking Contributions 

Under renormalization, there are continuum-hke mixings of the desired 4-quark lattice 
operators [38] , but there may in addition be perturbative mixing with operators of incorrect 
spin and taste. Further, there are discretization corrections to our identification of the spin 
and taste of the operators, which can arise from variations in the heavy-quark field over a 
spatial cube, as well as higher order terms in the position-space spin-taste formalism [9]. 
However, incorrect spin-taste operators generated by any of these causes cannot contribute 
to nonanalytic terms in the matrix element until NNLO. 

In order to see this, consider the standard power-counting in (rooted) staggered chiral 
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perturbation theory {rSxPT): 

p'^ ^m^a^ . (80) 

The 5° — mixing four-fermion operators, when translated into the chiral effective theory, 
are of 0{1) in the aforementioned power-counting scheme. Thus their LO (tree-level) contri- 
butions to the relevant matrix elements are of 0{1), and their NLO (one-loop) contributions 
are of C(p^). 

Perturbative mixing with wrong-taste operators can occur at one-loop order in as- In 
the chiral effective theory, these wrong-taste operators would thus enter with coefficients 
of 0{as/4:n). As was shown in Ref. [39], which considered the contribution of wrong-taste 
operators to neutral kaon mixing in rSXPT, as/4:TT is numerically of the same size as the 
taste-breaking factor a^a| on the a ~ 0.12 fm MILC Asqtad ensembles [40]. Thus the 
appropriate way to include the strong coupling constant in the rHMSXPT power-counting 
is: 

~ m ~ ~ as/4:Tc . (81) 

One-loop chiral diagrams involving the wrong-taste operators from perturbative mixing 
would therefore contribute to matrix elements only at NNLO, 0{as/4:TT p^), higher order 
than I am considering here. 

Wrong-taste operators may also occur because of the 0{a) corrections to the taste iden- 
tification of the operators and the interpolating fields, Eqs. (40) and (44), coming either 
from the variation of the heavy quark field over the hypercube (see Eq. (36)) or from 0{a) 
corrections to spin-taste identification of Eq. (38). Since the matrix elements are taste 
conserving at tree level in rHMSXPT, any such taste-violating effects must appear twice, 
inducing 0(0^) corrections. These effects can then be absorbed into a^-dependent analytic 
terms (see Eq. (95) below). Nontrivial terms could appear at one loop in chiral perturbation 
theory, but the extra factor of implies that such terms are again effectively NNLO. 

Consequently, only the LO taste breaking in the 4-quark operators, together with taste- 
breaking terms in the LO pion chiral Lagrangian, will modify the one-loop continuum chiral 
logarithms, and these modification are what have been calculated above. Note however, 
that in more highly improved versions of staggered quarks, taste violations (the terms in 
Eq. (81)) will be reduced relative to the Asqtad case, but perturbative mixing may not be 
similarly reduced. In such cases, it may be more reasonable to consider the 0{as) pertur- 
bative corrections to be of LO in the chiral expansion, and their one-loop chiral corrections 
to be NLO, the same order as the corrections computed in Sec. IV B. I make some more 
comments on this point in Sec. VI. 

V. FINAL RESULTS 

We can now combine results from different graphs to get the complete tadpole and sunset 
contributions. For the correct-spin tadpole contribution, adding Eqs. (68), (69), (72) and 
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(74) gives 



ie{Mf' } 



dixy 



(3)- 



i6{M(?'} 



+ (l^ ^ A) 



^2) 



dX 



xy 



Xy 



-E bSI({a^S};W})i 



(83) 



The incorrect-spin tadpole contributions come from Eqs. (73) and (75). Adding them, 
and using the values of [i^f^ in Table I and of ({k,,p) in Table II, we have 



- - 



5IxP - AXxA + 18Xx T - 4Xx,v - 5Xx,/ + 



V 



2 1-/ 
—a 



2(/32 + /33) 
/3i 



E^S({mS};{M)x 



E bKU{A/S};{^.})I 



, (84) 
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(87) 



A' 1 16 
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Xx P + 4Xx 4 - 6Xx,T + 4Xx V - Xxj + 
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X^,P-Xx,,-2a^<5;. 



ax 



Xv 



-E ^SU{^S};{M):^.,v 



- 2aH' 



(3)l. 



ax 



XA 



-E 



(3)i. 



In the continuum limit, when all taste propagators become degenerate and o^^'va ~^ 0> 
wrong-spin terms clearly vanish. 
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For the correct-spin sunset diagrams, adding Eqs. (76) and (78) gives 



(89) 



The incorrect-spin sunset contributions come from Eqs. (77) and (79). Again using the 
values of ({k, p) in Table II and of /^i"^ in Table I, we have 
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, (90) 



Q(2) 
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, (91) 
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Q(5) 



J TT 




njA* 
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xy 
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(3)- 



9'^X,A 



(94) 



In comparing Eqs. (62), (82), (83) and (89) to the continuum results of Ref. [11], one 
needs to be aware of the many differences in notation. In particular, the particles called X 
and TT in Ref. [11] are called r] and X, respectively, here. Taking the notational differences 
into account, it is straightforward to check that the continuum limits of the above equations 
reproduce the results in Ref. [11]. 
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The analytic terms in Eqs. (60) and (61) are of the form 



analytic terms = L^^^m^ + L^^\2mi + ruh) + hf^a^ , 



(95) 



where Lv and Ls are continuum low-energy constants, and La are lattice low-energy 
constants, summarizing the effects of taste-violating analytic chiral operators at NLO. As 
indicated, these constants depend on n, i.e., on the operator whose matrix element is being 
calculated. It is straightforward to see that these are in fact the analytic terms that may 
appear by considering the effects of adding mass or taste-violating spurious to the chiral 
operators in Eq. (57). The parameters L'f'\ and together with the parameters /3„ 
and /3'„ in Eqs. (60) and (61) are to be determined from the fits to lattice data. 

VI. CONCLUSIONS AND DISCUSSION 

I have calculated neutral B mixing to one-loop in staggered chiral perturbation theory 
for the complete set of Standard Model and beyond-the-Standard-Model operators, Eq. (1). 
My results are given by Eqs. (60) and (61), with expressions for the various terms listed in 
Eq. (62) and Eqs. (82) through (95). 

The construction of the operators on the lattice as local products of local heavy-light 
bilinears, coupled with the use of staggered light quarks, results in the appearance of "wrong 
spin/taste operators" that are 0{1) in the lattice spacing. Their contributions to the matrix 
elements considered here are suppressed to NLO because violation of taste symmetry is 
required. At that order, they induce mixing of the operators, summarized in the quantities 
T^"), Eqs. (84) through (88), and Q^'^\ Eqs. (90) through (94). However, as long as all five 
On are analyzed simultaneously, there are no new low-energy constants induced by these 
effects: the constants /3„ and (3'^ are all already present in the continuum. 

Effective operator mixings may come from three additional sources, weak-coupling pertur- 
bative corrections, corrections to the position-space spin-taste construction, and corrections 
to the taste identification of the operators and interpolating fields. I argue that any nonan- 
alytic terms that arise from such mixings are effectively NNLO, higher order than what has 
been considered here. 

Because the light staggered quark is converted to a naive quark in the lattice represen- 
tatives of the operators, the relationship between the staggered and naive quark fields plays 
an important role in my analysis. In particular, a naive quark is equivalent to four copies of 
staggered quarks, and the resulting "copy symmetry" can be used to simplify the calcula- 
tions. An interesting resulting feature is that the taste structure of the operators is not the 
same for all diagrams, but depends on the quark flow. 

Lattice computations in Refs. [1, 2] focused on the calculation of the quantity ^ = 



ral effects of the wrong spins were not known at the time of the HPQCD calculation [1] and 
were therefore omitted from the analysis and error estimate. In the Fermilab/MILC calcu- 
lation [2], the complete rHMSXPT expressions were available, but the matrix elements of 
operators other than Oi were not calculated, preventing a direct inclusion of the wrong-spin 
effects. However, it was possible to estimate the error of omitting these effects by using a 
small subset of new data to investigate the other matrix elements. The result, = 1.268(63) 




comes from the matrix element of operator Oi. The chi- 
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included a 3.2% error from this effect, wliich was the second largest source of error. In the 
ongoing second-generation Fermilab/MILC project [41], matrix elements of all five operators 
On are being computed, which means that the complete rHMSXPT expressions can be used 
in the analysis, and there will be no "wrong-spin error." Of course, a chiral/ continuum 
extrapolation error will remain. 

For future lattice computations of mixing with "highly improved staggered quark" (HISQ) 
ensembles [42], taste violations are sufficiently reduced that the power counting used here, 
Eq. (81), may no longer be appropriate. Depending on the range of lattice spacings studied 
and the statistical errors in the data, taste violations may in fact be so small that continuum 
heavy meson XPT might prove adequate for describing the data. More likely, one will want to 
use the rHMSXPT forms calculated here, but it may be necessary in addition to include the 
NLO chiral effects of the operators that enter through weak-coupling perturbative mixing, 
since such effects may no longer be much smaller than taste- violating NLO effects. Including 
such effects in NLO rHMSXPT will however be straightforward, since the chiral logarithms 
for the complete set of operators, Eq. (1), have already been calculated above. 
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Appendix A: B Parameters 

It is sometimes convenient to express the mixing matrix elements in terms of B (or "bag" 
parameters. A fairly common set of definitions is given for example in Ref. [11]: 



(bI\o:\bI) 



(1) 



r^nR'MlJl^B^^ forn = 2,3,4,5 



(Al) 
(A2) 



where Jb^ is the decay constant of the meson, Mb^ is its mass, R = MB^/{mh + m^), 
and ri2 = —5/3, rj^ = 1/3, 774 = 2, and 1]^ = 2/3. Relativistic normalization of the states is 
assumed in these expressions. 

The expression for the decay constant in rHMSXPT, including heavy-meson hyperfine 
and flavor splittings, is given in Eq. (6.20) of Ref. [7]. To convert it to the current notation, 
we just must replace the chiral logarithm functions J and i with T-t and X using Eqs. (63) 
and (64). Using quantities defined above in Eqs. (62), (68) and (69), we may write the result 
as 



1 + 7: 



'x.Fig.S "T 'x,Fig.4 



-(") 



xb 



(A3) 
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It is not surprising that W^^, T^^^pigS' ^^"^ Xl-mgi appear, because the wave function and the 
tadpole contributions of Figs. 3 and 4 are factorizable: they affect only one meson and one 
bilinear of the 4-quark operators, so they are exactly the contributions that appear in the 
decay constant. 

As always, heavy-light chiral perturbation theory is expressed in terms of the nonrela- 
tivistically normalized states of heavy quark effective theory. Thus the low energy constant, 
$0; that describes the decay constant at tree level in chiral perturbation theory includes a 
factor of -^/Ms^. Similarly, with the relativistic normalization used Eqs. (Al) and (A2), one 
factor of Mb^ needs to be included in our expressions for these matrix elements in terms 
of the parameters /3„ and /3^. Taking these normalization factors into account, and using 
Eqs. (60), (61) and (A3), we have 



B 
B 



(1) 

{"■) 



/3 



(8/3)^0 



^ (l + 5^ + f}^^ + Qi^) + Qi^)) + analytic terms 



(A4) 



(l ± 5. + ri")) + (si") + Qi")) + analytic terms , (A5) 



where n = 2, . . . , 5 in the second equation, and the upper (plus) sign is for n = 2, 3, while 
the lower (minus) sign is for n = 4,5. Here the wave function and tadpole contributions of 
Figs. 3 and 4 have canceled, and Sx is a new correct-spin tadpole term that comes only from 
the nonfactorizable tadpole diagrams. Figs. 5 and 6. Combining Eqs. (72) and (74) gives 



+- 



(2)i 



dl 



x,i 



(A6) 



The other chiral logarithm functions in Eqs. (A4) and (A5) are given in Eqs. (84) through 
(94). The analytic terms in Eqs. (A4) and (A5) have the same form as in Eq. (95) (with, of 
course, redefined low energy constants). 
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